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Biophysical  aspects  of  cyclodextrin  interaction 
with  paraoxon 

Sunil-Datta  Soni,^  Jayendra  B.  Bhonsle*^  and  Gregory  E.  Garcia^^ 


Cyclodextrins  are  torus-shaped  polymers  of  glucose  that  can  bind  organophosphorous  compounds  such  as  nerve  agents  and 
pesticides.  We  demonstrate  here  that  cyclodextrin  can  bind  up  to  two  paraoxon  molecules  with  a  of  6775  IVr\  Molecular  modeling 
shows  that  the  paraoxon  appears  to  bind  in  polar  opposite  orientation  and  have  an  average  binding  energy  of  —89  Kcals/mol. 
Published  2013.  This  article  is  a  U.S.  Government  work  and  is  in  the  public  domain  in  the  USA. 
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Introduction 

Because  of  the  initial  discovery  that  cyclodextrins  (CDs)  can  bind 
and  catalyze  the  hydrolysis  of  organophosphorus  (OP)  pesticides 
and  nerve  agents/^”^'  these  compounds  have  been  of  interest  as 
potential  scavengers  of  OPs.  The  intrinsic  ability  of  CDs  to  hydrolyze 
OPs  has  been  improved  by  several  modifications  such  as  capping 
on  the  primary  rim  and  addition  of  nucleophile  adducts  on  the 
secondary  rim.^'*'^'  To  make  improvements  in  CD  design  and  gain 
insights  into  the  structural  relationships  to  OP-scavenging  activity, 
we  began  by  examining  the  interactions  of  fl-CDs  with  paraoxon 
(PX)  by  NMR  techniques  and  molecular  modeling.  Whereas  NMR 
techniques  have  been  used  to  evaluate  the  interactions  of  CDs  with 
other  molecules,  as  far  as  we  know,  this  is  the  first  report  on  the 
interaction  of  CDs  with  this  pesticide  using  this  technique. 

Materials  and  Methods 

The  fi-CD  and  paraoxon  were  purchased  from  Sigma  Aldrich, 
Saint  Louis,  MO,  and  used  without  further  purification;  purity 
was  checked  using  proton  NMR  and  was  taken  into  account  in 
the  calculations.  The  ^-CD  was  98%  pure  as  supplied,  and  the 
water  content  was  15%.  All  compounds  were  dissolved  in 
99.8%  D2O.  The  NMR  experiments  were  performed  on  Varian's 
Unity  Inova  at  600  MHz  using  a  5-mm  Penta  PPG  probe  at  25  “C. 
Either  internal  4,4-Dimethyl-4-Silapentane-l -Sulfonic  acid  (DSS)  or 

a  b 


Figure  1.  Compound  structures,  a.  A"Cyclodextrin;  b.  paraoxon.  The 
number  schemes  are  shown  for  the  compounds. 


Table  1.  Chemical  shift  values  (ppm)  of/i-cyclodextrin  and  paraoxon 
at  25  °C 

Chemical  shift  (ppm) 


Compound  Resonance  Without  complexation  1 : 1  complex 


A-CD 

HI' 

4.938 

4.956 

H2' 

3.516 

3.536 

H3' 

3.834 

3.839 

H4' 

3.454 

3.475 

H5' 

3.733 

3.726 

H6',  6" 

3.750 

3.764 

Paraxon 

CH3  (triplet) 

1.361 

1.250 

CH2  (multi) 

4.332 

4.226 

H3,5  (doublet) 

7.445 

7.350 

H2,6  (doublet) 

8.334 

8.226 

residual  Deuterium  hydroxide  (HOD)  peak  was  used  as  proton 
reference.  A  standard  (s2pul)  pulse  sequence  was  used  and  90° 
pulse  calibrated  on  the  sample  with  a  3  s  relaxation  delay.  A  total 
of  32  scans  were  collected  for  each  data  set.  Data  were  processed 
using  MestReNova  software  (Mestrelab  Research,  Spain)  with  1  Hz 
line  broadening  and  zero  filling  to  32  K.  In  case  of  multiplets,  the 
average  chemical  shifts  in  parts  per  million  are  reported.  All 
proton  assignments  were  confirmed  using  2D  data  set  (COSY  and 
TOCSY). /?-CD  and  PX  were  mixed  and  equilibrated  to  25  °C  before 
proton  NMR  data  were  collected  and  chemical  shifts  were 
extracted.  The  structure  of  PX  was  built  and  minimized  with  MMFF 
force  field  and  charges.  The  molecular  volume  and  molecular  Van- 
der-Waals  surface  area  (measured  in  cubic  Angstroms  and  square 
Angstrom)  were  determined.  The  crystal  structure  of  J3-CD  (PDB 
ID=1DMB)  was  downloaded  (Protein  Data  bank  www.pdb.org) 
and  cleaned  as  follows.  The  fi-CD  ligand  was  extracted  from 
IDMB  and  minimized  with  MMFF  force  field  and  charges.  This 
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Figure  2.  NMR  analysis  of  paraoxon  (PX)  and ^-CD  interaction.  Job's  plot  analysis  (continuous  variation  method)  was  performed  for ^-CD  HI',  H2',  and 
H4'  protons  and  is  shown  in  a-c  respectively.  The  PX  CH3,  CH2,  H3,5,  and  H2,6  plots  are  shown  in  d.  (5  =  (5(complex)-^(free),  and  the  Mole  fraction 
r=  [CD]/([CD]  +  [OP])  where  CD=/i-CD  and  OP  =  PX.  The  Benesl-Hildebrand  plot  for  PX2-CD  Is  shown  In  e.  All  concentrations  are  In  millimeters. 


[Gy[H] 


[G]1H] 


Figure  3.  Plots  of  the  change  of  chemical  shift  association  constant  generated  from  titration  data  between  /9-cyclodextrin  and  paraoxon.  The  plots 
show  observed  and  calculated  data  for  the  ^-cyclodextrin  sugar  resonances  analyzed  using  nonlinear  regression  analysis  for  a.  HI',  b.  H2',  c.  H5', 
d.  H2  H8,  and  e.  H3  H5. 


wileyonlinelibrary.com/journal/mrc 
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was  used  for  the  docking  employing  the  FlexiDock  algorithm  on 
the  basis  of  genetic  algorithm  as  implemented  in  Tripos  Sybyl-X 
Version  1.3  molecular  modeling  software. 

Results  and  Discussion 

The  structures  of  /?-CD  and  PX  are  shown  in  Fig.  1.  The  chemical 
shifts  (ppm)  for  the  hosts  (/?-CD),  guests  (PX),  and  1 : 1  complex 
are  reported  in  Table  1.  The  NMR  PX/CD  data  were  compiled  to 
generate  values  to  the  Job's  plots  shown  in  Fig.  2.  The  continuous 
variation  method  [(FI)  +  (G)  =  (constant)]  of  Jobs's  plots  are  shown 
in  Fig.  2a-d.  The  AS  chemical  shifts  of  /?-CD  Fll',  FI2',  and  FI4'  are 
shown  in  Fig.  2a-c  and  those  of  PX,  CFI3,  CFI2,  FI2,5,  and  FI3,6 
are  shown  in  Fig.  2d.  The  fi-CD  plots  yield  maxima  of  0.33,  and 
the  PX  plots  yield  maxima  of  0.66,  indicating  that  more  than 
one  PX  molecule  is  binding  to  each  CD,  as  typically  a  1 : 1  ratio 
would  yield  a  maximum  of  0.5.  Additional  analysis  by  Benesi- 
Flildebrand  (B-FI)  plot  is  shown  in  Fig.  2e.  This  plot  did  not  yield 
a  straight  line,  further  indicating  that  multiple  PX  molecules  are 
binding  to  each  CD,  because  a  simple  1 : 1  complex  of  CD  and 
PX  would  yield  a  straight  line  from  which  the  association  constant 
Xav  would  be  the  slope.  Thus,  in  our  case,  the  association  constant 
was  determined  as  described  by  Flirose^®  where  FI  was  the  host 
(/?-CD)  and  G  was  the  guest  (PX).The  data  were  plotted  against 
[G]/[FI]  to  yield  the  plots  shown  in  Fig.  3.  The  calculated  average 
association  constants  were  derived  from  the  graphs  of  Fig.  3. 
The  observed  and  calculated  chemical  shift  values  from  nonlinear 
regression  analysis  were  subject  to  statistical  analysis  using 
SoIvStat'^'  and  are  reported  in  Table  2.  Because  of  severe  overlap 
of  some  of  the  fi-CD  proton  resonances,  only  some  of  the  values 
are  reported.  The  average  association  constant  (Xgv)  is  6557  M"' 
with  an  average  standard  deviation  of  1118  and  mean  of 
0.9932.  We  then  examined  fi-CD  accommodation  of  PX  by  molecular 
modeling.  The  computed  molecular  volume  and  molecular 
Van-der-Waals  surface  area  for  PX  were  found  to  be  779.56 
and  495.1 4  A^  respectively.  The  computational  modeling  of 
PX  into  the  /?-CD  cavity  is  shown  in  Fig.  4.  The  computational 
fitting  of  a  single  PX  into  the  fi-CD  cavity  is  depicted  in  Fig.  4a 
and  b  showing  the  wire  rendering  of  the  primary  face  view  and 


Table  2.  Results  from  the  nonlinear  regression  fit  of  HI',  H2',  H5',and 
H2  H8  and  H3  H5  shown  in  Fig.  3  as  association  constant  and  with 
statistical  errors 

Proton 

Parameters 

<5h 

<5GHG-complex 

^av  M 

HI' 

Values 

4.9359 

5.0072 

6659 

SDs 

0.0005 

0.0023 

675 

R",  SE(y) 

0.9981 

0.001 

— 

H2' 

Values 

3.5146 

3.5873 

7866 

SDs 

0.001 1 

0.0033 

1384 

R",  SE(y) 

0.9930 

0.0023 

— 

H5' 

Values 

3.4528 

3.5278 

9276 

SDs 

0.0012 

0.0031 

1592 

R",  SE(y) 

0.9930 

0.0024 

— 

H3,  H5 

Values 

8.2038 

8.2885 

4971 

SDs 

0.0013 

0.0062 

1117 

R",  SE(y) 

0.9905 

0.0028 

— 

H2,  H8 

Values 

7.3385 

7.3916 

5013 

SDs 

0.0008 

0.0035 

1029 

R",  SE(y) 

0.9921 

0.0016 

— 

Mean  Kav 

6757 

a  b 


Figure  4.  Molecular  models  of  single  and  double  paraoxons  (PXs)  within 
fi-CD.  Wire  figures  of  single  PX  binding  are  shown  in  a  (primary  face)  and 
b  (secondary  face).  The  wire  figures  of  double  PX  binding  are  shown  in  c 
(primary  face)  and  d  (secondary  face).  The  lateral  views  of  single  e  and 
double  f  PX  binding  are  shown  with  the  primary  face  on  the  left  and  PX 
rendered  as  capped  sticks.  Oxygen  is  rendered  in  red,  phosphorus  in 
orange,  and  nitrogen  in  blue. 

secondary  face  views.  The  docking  of  two  PX  molecules  into  the 
fi-CD  cavity  in  the  wire  rendering  form  is  shown  in  Fig.  4c  and  d 
for  the  primary  and  secondary  face  views.  The  docking  energy  of 
the  two  PX  in  the  PX2/CD  is  — 89Kcal/mol,  while  for  one  PX  in 
PX/CD,  it  is  —55  Kcal/mol.  In  the  lateral  face  views  shown  in  Fig.  4e 
and  f,  the  PX  and  2PX  molecules  are  shown  as  Stick  renderings  to 
highlight  them  vis-a-vis  the  fi-CD  molecules.  The  model  shows  that 
the  phosphate  and  nitro  portions  of  the  PX  molecule  protrude  out 
of  the  fi-CD  cavity.  For  2PX  binding,  the  phosphate  and  nitro  of  the 
second  PX  are  polar  opposite  the  other  PX.  This  arrangement  would 
allow  the  aligning  of  the  phenyl  portions  that  may  produce  tt-tz 
interactions  that  enhance  the  docking  energy  for  the  second  PX, 
although  this  would  not  explain  the  large  difference  in  docking 
energies  as  yet  other  unidentified  interactions  could  contribute. 
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